CHEP-PKU/1/01-2008, IMSC/2008/01/01 



Spontaneous symmetry breakdown in fuzzy spheres 

C.R. DatQ 

Center for High-Energy Physics, Peking University, Beijing 100 871, China 

S. DigaH and T.R. Govindarajaifl 

Institute of Mathematical Sciences, C.I.T. Campus, Taramani, Chennai 600 113, India 

We study and analyse the questions regarding breakdown of global symmetry on noncommutative 
sphere. We demonstrate this by considering a complex scalar field on a fuzzy sphere and isolating 
Goldstone modes. We discuss the role of nonlocal interactions present in these through geometrical 
considerations. 

PACS numbers: 12.60.Rc; 12.10.-g; 14.80.Hv; 11.25. Wx; 11. 10. Hi 



I. INTRODUCTION 

Recently there has been interest in the non-perturbative numerical studies of quantum fields on 
fuzzy spaces. This was motivated by the fact that implementation of noncommutative geometry 
appears through novel features in such models. For example, the IR/UV mixing, even though 
absent in the finite matrix setting appears through an anomaly which generates it in the continuum 
limit. This suggests the possibility that fuzzy spaces may serve as novel regulators of field theories 
and bring out new features which were absent in the conventional lattice regularisations. This 
program falls well within the front line of research activity of study of noncommutative geometry 

iiiiiaiitatni. 

Earlier studies on simulations have focused on scalar fields on fuzzy spheres ( S%) and led to 
the demonstration of new phases characterised as nonuniform ordered phase @, H (Mini. In the 
continuum infinite volume limit this nonuniform phase appears as stripe phase in the Groenewold- 
Moyal space [T3 . [l3| . These nonuniform phases are related to the existence of meta stable states and 
we have demonstrated this in our earlier analysis using what is known as pseudo heat bath method 
[TT[ | . We could also show the existence of meta stable states as well as their connections to many 
phases. 

Those new "stripe" phases were conjectured by Gubser and Sondhi [12], where they pointed 
out that this translation non-invariant phase would exist only in dimensions d > 3. However the 
numerical simulation of Ambjorn and Catterall show the existence of such a phase even in d = 2 
This seems to be contradicting Coleman-Mermin- Wagner (CMW) theorem that continuous 
symmetries cannot be broken spontaneously in d = 2 [l4j. But CMW theorem can be expected 
to be valid only in local field theories [l5| . On the other hand NC field theories are inherently 
nonlocal and can be expected to show new features. However Gubser and Sondhi while analysing 
NC field theories point out that the infrared problems are worse in these theories and hence CMW 
theorem cannot be violated. But numerical simulations show the results to be other way. In this 
connection it is important to see whether long range order or symmetry breaking can be seen for 
global symmetries other than space time symmetries. This assumes importance, since Gubser and 
Sondhi explicitly use O(N) and 0(2) symmetries in showing their results. 
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We analyse in this paper this important question, namely what happens to CMW theorem and 
long range order in the fuzzy spaces 15]. We consider a complex scalar field theory with a global 
U(l) and study the full implications of the noncommutative nature of the underlying geometry. As 
pointed out earlier spontaneous symmetry breakdown (SSB) and long range order are obstructed 
in 2D theories by the infrared divergences associated with Goldstonc modes. While nonlocal field 
theories can escape these conditions, it is not obvious whether they do so in these circumstances. 

This paper is organised as follows: In Sec. 2 we describe the model, notations and address the 
questions. In Sec. 3 we look at the aspects of simulations and demonstrate the spontaneous breaking 
of U(l) symmetry thereby evading CMW theorem. In Sec. 4, we analyse questions of Goldstone 
modes, the role of non-locality and subtle issues in isolating effects of continuum limit as well as 
infinite volume limit. In the last Sec. 5, we present our results and conclusions, taking into account 
already existing theoretical studies. 



II. SSB ON 2D NC SPACE AND CMW THEOREM 



The standard action used for studying the U(l) in two dimensional (2D) Moyal spaces is given as 

S = j d 2 x [|<9$| 2 + r|$| 2 + Ai$*$*$**$* + A 2 $*$**$*$*] (1) 

Finite temperature behaviour of the $ field with fluctuations is studied by varying the mass param- 
eter r. One would expect the average of <& to take the form, 

in the mean-field theory. In commutative two dimensional space the U(l) Goldstone mode destroys 
the above condensate. So, there will not be SSB of the U(l) symmetry in these spaces. As discussed 
above, the status of U(l) SSB is not clear in 2D NC spaces. Very little has been done from the 
non-perturbative side to study this issue. 

We will present in this note that the spontaneous breaking of internal symmetries can be expected 
in Moyal space times when space-time symmetries are broken. Conventionally the obstruction to 
SSB in 2D will come through the infrared divergence from the resulting massless Goldstone boson. 
Gubser and Sondhi argued that the infrared behaviour is worse in 2D NC space and hence will 
continue to obey CMW theorem. But the simulations point out that this to be not the case. We 
will argue how this can be understood with specific requirements in these theories. 

To obtain SSB of U(l) we should have the complex field <E> in the ground state such that 

J $ d 2 x = and $ ^ . (3) 

This will translate to Tr<I> = in Moyal space. This implies that when $ is expanded in Fourier 
modes $(fc), the zero mode should be absent. That is, 

$(fc = 0) = 0. 

Correspondingly, on the sphere $ should not have a singlet component in the angular momentum 
basis. This naturally gives a cutoff k c fixed by the parameters of the effective action including 
quantum fluctuations. This cut-off in the integral avoids infrared singularities. The important 
factor responsible for this is the existence of states which are not translationally invariant, such as 
the stripe phase (or nonuniform phase). Hence U(l) can be broken spontaneously only along with 
translation symmetry violating nonuniform phase. We will see how these features appear in our 
simulation studies. 
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III. SIMULATIONS ON FUZZY S 2 



In this work we study this problem non-perturbatively on a fuzzy sphere. The fuzzy spheres are 
described by the coordinates obeying the algebra 

[x i ,x j ]= 7;^, , j: x ? = r2 - ( 4 ) 

v JU + 1) i 



Here i? is the radius of the fuzzy sphere, a / yj(j~- PT) measures the non commutativity in the fuzzy 
sphere. Here a is a function of R. On the fuzzy sphere the above action (Eq.l) reduces to, 



47T 

<?($) = -Tr 



^ | $] | 2 + R 2 (r|$| 2 + Ai^^^^t + As^^t) 



(5) 



Here the real and imaginary parts of <f> G Mat^r are N x N hermitian matrices. The quartic 
terms represent the self interaction of the $ field. These terms are crucial for the presence of 
noncommutativity and nonlocality in this problem. In the following we describe our strategy of 
studying the SSB of U(l) and the numerical simulations. 

The simulation of the above model for the SSB study of U(l) in fuzzy S 2 involves generating 
statistically relevant matrix configurations for a particular choice of parameters. The configurations 
are generated using "pseudo-heatbath" method, for details see 11]. The auto-correlation for larger 
TV simulations turn out to be large. To reduce the autocorrelation problem we use over-relaxation 
algorithm. However the number of over-relaxation steps required grows with N rapidly. This makes 
it difficult to get any reasonable statistics for large N. 

The CMW theorem and SSB for finite systems is subtle. For a finite system there will be tunnelling 
between all possible vacua, though the tunnelling rate is exponentially suppressed as the size of the 
system increases. Because of the tunnelling, the state of the system respects the symmetry. However 
this state of the system is unlike the symmetry restored state at high temperatures. In this situation 
ensemble average of the magnitude of the "magnetisation" is usually taken to be the condensate 
(order parameter) . Using this prescription we look for SSB of continuous symmetry for finite systems 
in 2D. In order to numerically prove the CMW theorem one has to show that the condensate vanishes 
in the infinite volume limit. The numerical study of finite size effects in commutative 2D show that 
the condensate vanishes logarithmically. One should expect that there will also be finite size effects 
for noncommutative spaces. So the finite R, N effects must be separated out to settle the issue of 
U(l) SSB on fuzzy sphere. There are various limits possible consistent with (R, N) — > oo. We do not 
consider the limit (i?, N) — > oo in which the resultant space becomes commutative and there should 
be no SSB of U(l). When the ratio R 2 /N remains finite for (i?, N) — > oo we have a thermodynamic 
limit which still maintains the non-commutative character. 

The numerical simulations always consider a finite N. As discussed above a non-zero (<E>) will not 
validate or invalidate the CMW theorem. We need to study the behaviour of (<E>) for various TV to see 
SSB of U(l) or otherwise. While we vary N the ratio R 2 /N is kept fixed. To simplify the numerical 
work we choose Ai = A2 = A/2 from the parameter space. The simulations were then done for a 
fixed A,r. The observables we measure are action S, and fluctuations of various elements of $. 
Each element &ij now is a complex number, so, 



'I'-.- >!,,<"' • (6) 

Note that because of the U (1) symmetry $^ will have a circular distribution in the complex plane, 
with (r)ij) being the average radius of the distribution. A non- vanishing (r)ij) in the N — > 00 limit for 
at least one element of $ will clearly signal SSB of U(l). In the following we discuss the simulation 
parameters, results and discussions. 
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IV. RESULTS AND DISCUSSIONS 



The parameter values we choose for our simulations are A = 1.0, r = —8. We have also repeated 
our simulations for r — —16. For these choices of parameters various values of N starting from 
N = 16 onwards are considered. The run time to get reasonable statistics grows rapidly with N. At 
present the largest value of N we have simulated is N — 100. We also have a few set of data with 
less statistics for N > 100. 

The simulation results showed that the off-diagonal elements of $ are always fluctuating around 
zero. The change in the behaviour of these elements with N is not so prominent. So, in this paper we 
show only the behaviour of diagonal elements of In Fig.l, we show the distributions of $n, $22 
and <i> mm , where mm correspond to the middle diagonal element or one of the middle diagonals if 
N is even. 




FIG. 1: 



From the figure we can clearly see that the distribution has a U(l) symmetry. The distributions 
of pair of elements such as and $jv-i+i,iv-i-i-i are found to be the same as expected. In order 
to show finite size effects we consider N = 19 and N — 64 in the above figure. We see that for all 
N, the radii of the distributions decrease as we move towards the middle elements. The difference 
here is that for N — 19 the radii are close to each other, but for N = 64 they are different by large 
amount. Indeed for N = 64 the middle diagonal distribution has a peak around zero in the complex 
plane. The difference in the property of <E> between N = 19 and N = 64 is purely finite volume 
effect. Also due to finite volume effects even the radius of distribution of &h($nn) has increased 
from N = 64 to N = 19. 

One of the main difference between commutative and non-commutative 2D space is the presence 
of degenerate/meta-stable states in the latter. In the numerical simulations one always chooses 
an initial configuration. For commutative 2D, whatever may be the initial configuration, there is 
always a unique thermalized state. The situation is completely different in the case of fuzzy sphere. 
Depending on the initial configuration we get different thermalized states. 

The detailed analysis of $ reveals that the phase of all diagonal elements is perfectly correlated 
for smaller N. For example in Fig. 2 we plot the phase of $n vs phase of §nn for N = 19. 

We find that each configuration for N — 19 the phases of all diagonal elements are fluctuating 
around same average value. Since the different radii for the distributions are not much different the 
form of for N = 19 is close to the form $ ~ le t6 . In Fig. 2b we plot the phase of $n vs phase of 
$atat for N = 64. Here, the correlation between 9n and 6nn is not as strong. These results were 
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(a) Distribution of 9\\ vs fljvjv for N = 19 and po ~ !• 




1 2 3 4 5 6 



©11 

(b) Distribution of flu vs Sjvjv for A" = 64 and po ~ 1- 



FIG. 2: 



obtained with a initial configuration <I>o = pol, with po ~ 1> With the initial condition, po ~ 0, we 
get same results for A = 19 as before, but for N — 64 the results are dramatically different. In Fig. 3 
we show the distribution of $ and 9%i vs #nn correlation for AT = 64. It can be seen from Fig. 3 the 
correlation between On and Onn is stronger in this case. 
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FIG. 3: Distribution of 9\x vs 6nn for N = 64 and po ~ 0. 




The difference in results for A?~ = 64 shown in Fig. 2b and Fig. 3 are coming from two (meta)stable 
states available for the system. Such (meta)stable states are absent for the same theory in commu- 
tative 2D space. For all different initial configurations for A^ > 25, 9n and 8nn were never in phase 
like Fig. 2a. As a result Tr($) decreases for larger N. For larger N average of Tr($) is found to be 
zero. This is expected as only these states with Tr($) = can survive the thermodynamic limit. 

From N = 19 to N = 64 there is transition in the form of $. Also the average of 7711 is smaller 
for N = 64. This result clearly shows that uniform ordered phase is not a stable state in the limit 
N — > 00. The transition in the form of 4>, when N is increased takes place around N — 25 for the 
choice of A, r in our simulations. This can be seen in Fig. 4a. 

The above results show that (rjn) seems to saturate for N > 40. In Fig. 4b, we show the histogram 
distribution of f?n for N = 40,64 and 100. The peak of the distribution hardly changes with N. 
On the other hand, the width for A^ = 100 is clearly smaller. In other words, the fluctuations of 7711 
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(a) N dependence of 1711,7722 and r\ mrn . 
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(b) Histogram of 7711 for N = 40, 64 and N = 100. 



FIG. 4: 



is sharply peaked. This suggests that the results for larger N will not deviate from what we have 
already seen. The results shown in Fig. 4a have been obtained with po ~ 1 for the initial configuration. 
We find that the results for 7711 or tjnn hardly change with different initial configurations. Even if 
there is a small change, their values do not show any N dependence for larger N. We caution, here 
that the statistics for N > 100 are very small and the autocorrelation is substantial. 

The saturation of (7711) beyond N — 40 up to the largest simulated N suggests that <£> does not 
vanish in the non-commutative thermodynamic limit. The form of $ we have obtained from the 
simulations also show that the ground state is not a uniform ordered phase, but rather a non-uniform 
ordered phase. This shows that all the generators, L x ,L y and L z are spontaneously broken in the 
ground state. This implies that there will be three Goldstone modes. 



V. CONCLUSION 



Having presented various simulations on fuzzy sphere with the intention to understand the issue 
of CMW theorem we will now summarise the conclusions. While considering QFT's in Moyal space 
time Gubser and Sondhi [l2| concluded the possibility of translation symmetry violating stripe phase 
in dimensions d > 3. Ambjorn and Catterall [l3[ have demonstrated this phase in d = 2 through 
Monte Carlo simulations. This was followed by studies on Sp and Sp x R [l(| EH [13 • Gubser and 
Sondhi also conjectured that internal symmetries like O(N) cannot be broken, as expected from 
CMW theorem. But CMW theorem explicitly uses locality of interactions and NC field theories 
are inherently nonlocal. We have argued that the Fourier transform of the complex field should 
be vanishing at zero momentum, as the important criterion for global symmetry breaking phases. 
This translates to vanishing of the singlet component of the field in the angular momentum basis in 
fuzzy spheres. In our simulations we find evidence for the breaking of global symmetries in 2D NC 
field theories. The existence of nonuniform stable states contribute at higher temperature, thereby 
breaking the translation and internal symmetries. Translation symmetry breaking stripe phase is 
required for obtaining the internal symmetry violation and long range order and these two occur 
together. 
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